AD-785  231 


A  UNIFIED  THEORY  OF  SOLID  PROPELLANT 
IGNITION.  PART  1.  DEVELOPMENT  OF 
MATHEMATICAL  MODEL 

H.  H.  Bradley.  Jr. 

Naval  Weapons  Center 
China  Lake.  California 

August  1974 


DISTRIBUTED  BY: 


fbtioul  Ttcttcai  Minutm  Stnrict 
(L  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Read,  SpriaffeW  Vs.  22151 


UNCLASSIFIED 


MCu*<Tv  0.*tK»lC*Ti&»  O*  T».t  »•«(  HWi  Dot  lawW) 


M2*££2L 


REPORT  DOCUMENTATION  PAGE 


V  4|^cAf  - 

wc  TF  5618  •  PART  1 


*eax>  wrrtucnoM 

BtfOM  COUPLE  TP*  G  fOgi 


1  SOVT  ACCCU40M  M> 


*  tiT^t  ‘m+  Mwrt*) 

A  UNIFIED  THEORY  OF  SOLID  FROPOUUfr  IGNITION- 
FAX!  1 .  DEVELOPMENT  OF  KATKDUTICAL  HOOEL 


A  Tl*C  0*  *  ACMOO  COVCRCO 

Development  Program 


H.  H.  Bradley,  Jr. 


4  oso.  report  «v«»CR 

rrsjmn  ©TEnsr  «usnvr7 


*  a*og*a-  ITcmcvT  project.  TI5T 

*M*  4  •cane  umiT  «UNMM 

01D  331-001/200-1/UR 
024-02-02 


t  0«64«|lX  ATlO*  •••€  4*0  AOO«(tt 

Naval  Weapon*  Center 
China  Lake,  Calif.  93553 


««  COiT»3LLi»,%  0»»iCt  ham  a«o  *OD«(U 


»*  REPORT  0At« 

AUCUST  1974 


~ 


*4  woo  Tokn»6  A<UaCv  »n»I  4  AOOACl^N  mttrnrnm §  •—  C 


‘1  «CU«*1T»  CWAU  f<4  AN  M»i*> 

UNCLASSIFIED 

II*  ^Mv^AT,oil>o^ckAO^# 


W  0»$**iO\jT.o«n 


•  T  rmt  mmi  A 


Approved  for  Public  Release,  distribution  unlimited 


»T  0«lYANA^Tio«i  ITATf  MCMT  f»  mm  inwrf  «•  0l«c*  J4.  If  AlAww  Irw 


*4  »UA#WCJI* 


>4  *(V  vO^Ot  fCmiA»»  «•  ft— f—  »Mt  if  ■unify  aN  mM%  l 

Ignition  Theory,  Solid  Propellant 


>  AStTMACT  fCmm 

See  back  of  fora 


w  ft—rw  AA  If  —  CMitry  •—  lAw-Ify  *r  Mm*  « 


"^NATIONAL  TECHNICAL 
INFORMATION  StRVICE 

V  S  Co 

W**im*l  taA  :im 


DO  1473  to.T.o.0*  »«>¥*»  «oMOc«Tt  •  UDCLASSIFIH) 

t/m  MiHifiwi  |. 


iicvStt  nSSncIrioS  o»  Twit  »A4«  rC C Kt# 


WCLASSOTP _ 

cimryicitioii  Of  tun  cm  So*— w 


(V)  A  Unified  Theory  of  Solid  Pro¬ 
pellant  Ignition.  Part  I.  Development 
of  Mathematical  Model ,  by  H.  R.  Bradley, 
Jr.  Chlnal  take,  Calif.  Ravel  Weapons 
Center,  August  1974  ,  36  pp.  (KWC  TP 
5618,  Part  1,  publication  UNCLASSIFIED.) 

(U)  The  governing  equations  of  solid 
propellant  Ignition  Is  one  space  dimension 
are  derived  from  conservation  and  trans¬ 
port  considerations.  The  resulting  math¬ 
ematical  model  contains  as  subcases  the 
existing  gas  pt.ise,  solid  phase,  and  sur¬ 
face  Ignition  theories.  Hors  Important, 

It  allove  the  effect  of  the  site  of  chem¬ 
ical  reaction  on  the  Ignition  transient 
to  be  examined  with  consistent  and  compat¬ 
ible  boundary  conditions.  Surface  re¬ 
gression  and  the  first  order  effect  of 
oomanlty  Lewis  number  are  retained  to 
allow  examination  of  the  effect  of  these 
phenomena  which  are  often  neglected. 
Several  second  order  terms  are  dropped 
from  the  more  general  equations  In  the 
Interest  of  mathematical  traceability; 
these  Include  variable  specific  heat  and 
transport  coefficients.  Finally,  the 
physical  equations  are  transformed  for 
computer  coding.  The  computer  program 
and  results  of  a  parameter*  stud*  will  be 
published  as  separate  reports. 


JL _ 

UNCLASSIFIED 


MCv*«TV  ClAM4PtCATiO«l  OF  TMIft  PUCTMm  Dm  Sw—  4} 


1 


n,  i-r-y 


it&'i.S. 


mm 


•a,  ifowaraj  and«r  Iml  Ordanc*  Syr 
on  3JI-O0I/20O-1/U*  OT4-O2-02.  A  yr 
rt  «M7.  wens  rr  »J4),  —Tim  tte 


Fnbliafaad  by  . 
Collation  .  . 
first  printing 


irch  Department 
.  .  .  Corn,  b  Imtm 
.175  numbered  copies 


1*1 

/II* 


KJC  IT  5618,  ?«rt  1 


COWTDtTS 

tatocUturt .  2 

1.  Introduction  .  6 

2*  Physical  Description  of  Problea  •  •  . . .  •  8 

3.  Mathematical  Formulation  ........  .  .  •  .  10 

4.  Manipulation  of  Equations . 26 

Appendixes: 

A.  Interface  Conditions  for  Equilibrium  Pyrolysis  of  One 
or  More  Solid  Phase  Constituents  .  31 

I.  Simplified  Analysis  of  Effect  of  Absorption  of  Had  lent 
Ensrpy  by  Css  Phase  Species  •  •  •  . . 33 


WC  V  5618.  f«rt  1 


MKENCLAttftl 

Af  Constants  In  Eq.  79.  Saa  Kay.  following  this  Noaanclatura 

for  defining  expressions 

*i'*i  Stoichiometric  coefficients  for  bulk  and  surface  reactions 
Defined  in  Section  3*2.2 

CsCjy  Specific  heat  and  average  specific  heat  per  unit  ease 
c  Molar  concentration 
0,V  Coefficients  of  mass  diffusion 
d  Total  de*rlvatlve  operator 

fj  Activation  energy.  Explanation  of  numerical  subscripts  given 
in  Section  3.2*2 

F7J$  Reaction  orders  of  fuel  in  gas  phase  reactions 

M  Total  specific  enthalpy  per  unit  mass,  including  chemical  and 
thermal 

h  Thermal  specific  enthalpy  per  unit  mass 
bf  Reference  enthalpy 
Im  Unattenuated  initiating  flux  density 
j  Diffusion  mass  flux  density 

Jtj  Mass  source.  Explanation  of  numerical  subscripts  given  in 
Section  3.2.2 

tie  Lewis  number  ■  p CV/X 

Surface  mass  flux  density.  Explanation  of  numerical  subscripts 
given  in  Section  3.2.2 

07,08  Reaction  orders  of  oxidiaer  in  gas  phase  reactions 
p  Pressure 
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Ol  Ruts  of  reaction,  positive  0  being  exothermic.  Subscripts  1, 

7  and  8  correspond  to  ructions  1,  7,  and  8  (see  Sections  3.2.2.1 
end  3.2.2. 2).  Subscripts  5,  6,  10,  11,  and  12  correspond  to  the 
•quivalent  subscripts  on 

q  Total  thermal  flux  density  relative  to  use  average  phase  velocity 
qr  Radiant  hut  flux  density 

qM  Surface  rate  of  hut  evolution 

R  Universal  gas  constant 

rj  Chemical  reaction  rate  (see  Sections  3.2.2. 1  and  3. 2.2.2) 

Sj  Surface  sources .  See  Key  for  defining  expressions 
55*56  Reaction  orders  of  oxidizer  in  surface  reactions 
T  Temperature 
t  Time 

V  Total  specific  internal  energy  per  unit  mass 
u  Velocity  relative  to  stationary  coordinate  system 
Vi  Volumetric  sources.  See  Key  for  defining  expressions 
v  Maas  average  velocity 
Y£  Dlffuslonal  velocity  of  species  i 
Xj  Hole  fraction  of  species  i 

x  Space  coordinate 

t£  Mass  fraction  of  species  1  (aee  Eq.  42  and  40-54  for  subscript 
munlngs) 

Yfc  Initial  solid  fuel  mass  fraction 

yo*  Initial  gaseous  oxidizer  mass  fraction 

Zj  Generalized  dependent  variable;  Includes  mass  fraction  and 
enthalpy.  Zj  -  Tj  for  J  •  1,  2,  ...  12,  Z13  -  hQt  Zn  m  hc 
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o  Coordinate  transformation  parameter  (see  Eq.  73  and  76) 

0  Flux  attenuation  coefficient 
X  Coefficient  of  thermal  conductivity 
Uj  Molecular  weight  of  species  i 

vj  Pre-exponential  factor  in  rate  expressions.  Subscript  i  corresponds 
to  si&scrlpt  i  on  E 

o  Mass  density 

♦  Transformed  coordinate 

♦  Hovarth  transformed  space  coordinate 
7  Gradient  operator 

i  Partial  derivative  operator 
Subscripts 

c  Condensed  (solid  phase) 

?  Gas  phase 
s  Surface  (interface) 
o  Initial 


i  Referring  to  species  i  unless  otherwise  specified 
Key  to  numerical  subscripts 

»i  -  V'pi70  *  i  •**  9aa,,iJ 

«  z  i  m  20  f  lit  12  f  74 

Bf  -  Xc,  'c9  *  1  *  2'  "*9 

-  0  :  i  •  10.  U,  12 

-  Bi/Lo  :  i  -  13 

-  ac^c/Ce  *  1  *  14 

*»©  ilml,  3,7,  8,9 
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l.  umtooucnos 


Combustible  or  txploslvt  Ingredients  are  sn  essential  component  of 
ordnance  Items  ranging  in  t  «  from  small  arms  to  long  range  missiles* * 
Solid  propellant  charges  for  these  devices  consist  of  premixed  fuel  and 
ok  Miser  combinations  uhleh  must  remain  In  a  metastable  state  throughout 
the  useful  storage  life  of  the  Item*  Ignition  Is  the  stimulation  of  the 
meter la 1  by  any  of  various  energy  sources*  overcoming  the  metastable 
barrier.  The  term  Ignition*  as  used  In  this  report*  Is  primarily  con¬ 
cerned  with  such  stimulation  at  and  near  the  propellant  surface*  De¬ 
pending  upon  the  application*  there  may  or  may  not  be  an  Identifiable 
period  of  steady  state  operation* 

In  a  rocket  motor*  between  Initiation  and  final  burnout*  there  occurs 
*  series  of  Interrelated  processes  Including  mass  diffusion*  chemical 
reactions*  fluid  dynamics*  sod  beat  transfer  by  various  modes*  In  order 
to  simplify  the  Ignition  problem  end  to  place  the  current  Investigation 
In  proper  perspective*  It  is  Instructive  to  consider  three  phases  of 
Ignition*  First*  energy  from  the  initiator  Is  directed  toward  the  pro¬ 
pellant  surface*  The  ensuing  physicochemical  processes  lead  to  self- 
accelerating  reactions  with  evolution  of  gaseous  products*  Second*  the 
locally  Ignited  state  of  the  propellant  la  propagated  over  the  surface 
of  the  propellant  by  thermal  and  fluid  dynamic  means.  Third*  and  often 
simultaneously,  the  gases  evolved  from  the  Igniter  and  propel laut  produce 
a  pressure  build  up  In  the  device*  leading  to  some  desired  operation. 
Various  aspects  of  the  latter  two  phases  hsve  been  treated  elsewhere*** 


ds  Soto*  S.  and  H.  A** Friedman,  "Flame  Spreading  and  Ignition  Tran¬ 
sients  In  Solid  Grain  Propellants*"  AMEX  INST  AZJtOK  ASTRONAUT  J*  Vol*  3* 

Ho.  3  (March  1965),  pp*  405-12. 

*  United  Technology  Center.  Studies  on  Tg'.ition  and  Flans  Propagation 
of  Solid  Propellants ,  Final  Report  Covering  the  Period  22  June  1966  through 
22  June  1967 *  by  J.  D.  Kllgroe*  13  November  1967.  (UTC  2229-FX  prepared 
for  Jet  Propulsion  Laboratory*  NASA*  Pasadena*  Calif*) 

^  Peretx*  A.*  and  others. "Starting  Transient  of  Solid-Propellan*  Rockat 
Motors  with  High  Internal  Gas  Velocities,"  AMER  INST  AERON  ASTRONAUT  J  * 

Vol.  11,  No.  12  (December  1973),  pp.  1719-27. 

» 

*  Bradley*  H.  H.*  Jr.  "Theory  of  a  Homogeneous  Model  of  Rocket  Motor 
Ignition*"  (January  1964,  AIAA  Preprint  64-127),  6  pp. 
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and  are  not  considered  to  be  within  the  scope  of  the  current  study.  It 
Is  toward  the  theorttlcsl  septets  of  the  fir*t  phase,  which  Is  celled 
propellent  Ignition,  to  distinguish  it  from  rocket  motor  Ignition*  thet 
we  now  direct  our  ettentlon. 


1.1  BACXCROCKD  AND  REVIEW 

An  extensive  review  of  solid  propellent  Ignition  theories  wes  pub¬ 
lished  in  1966.5  It  was  noted  there  that  the  distinguishing  features  of 
the  theories  were:  (1)  the  assumed  site  (solid  phase*  gas  phase,  or 
Interfere)  of  the  key  chemical  process  leading  to  an  approach  to  steady-* 
state  combustion  through  a  self-accelerating  (due  to  Increasing  temper¬ 
ature)  and  self -Halting  (due  to  reactant  depletion  and  dlffuslonal 
effects)  reaction  rate;  (2)  the  nature  of  the  ignition  stimulus  (radiant 
•oergy,  hot  gas  or  solid,  or  hypergollc  chemical  reaction);  and  (3) 
choice  of  an  Ignition  criterion  (e.g.,  critical  temperature,  go/no-go, 
light  saissloo,  or  rate  of  taperature  rise). 

The  earliest  theory  chronologically  sod  the  simplest  mathematically 
was  the  solid  phase  thermal  theory  irfilch  was  an  outgrowth  of  thermal  ex¬ 
plosion  theory.  In  It,  a  solid  phase  homogeneous  reaction  Is  predicated 
(lustlfled  bv  known  solid  phase  exothermic  reactions  lu  double-base  pro¬ 
pellants)  with  the  Ignition  stimulus  consisting  of  s  radiant  source*  hot 
gas*  or  hot  solid.  An  obvious  deficiency  of  the  solid  phase  theory  was 
its  ability  to  predict  or  explain  fully  the  experimentally  observed 
effects  of  ambient  gas  conditions  (pressure  an d  composition)  upon  Ig¬ 
nition  time.  The  gas  phase  theory  wes  developed  In  connection  with  an 
experimental  program  In  which  a  solid  frel  was  ignited  by  high  tempera¬ 
ture  oxidizer  ^3  heated  in  a  shock  tube.  Several  versions  of  the  gas 
phase  theory  nave  been  presented  with  different  degrees  of  sophistication; 
all  contain  the  feature  that  fuel  vapors,  produced  by  heat  transfer 
from  the  hot  oxidizer  to  the  solid,  diffuse  into  the  gas  phase  and  sub¬ 
sequently  react.  The  Interface  theory  grew  out  of  experimental  obser¬ 
vations  in  which  a  highly  energetic  oxidizer,  viz.,  fluorine.  Ignited 
a  solid  fuel  hypergollcilly  at  room  temperature.  The  chemical  reactivity 
was  asstssed  to  be  confined  to  the  surface;  hence,  the  theory  has  been  re¬ 
ferred  to  as  bypevgollc  or  heterogeneous.  This  latter  terminology  should 
not  be  construed  as  denoting  relevance  to  heterogeneous  (composite)  pro¬ 
pellants,  although  results  of  the  analysis  have  been  extended  heurlstl- 
cally  to  Include  such  an  r  plication. 

In  addition  to  differences  in  formulation  of  the  three  theories, 
there  has  been  a  diversity  of  choice,  both  in  experiment  and  theory,  re¬ 
garding  the  attainment  of  Ignition.  Experimentally,  use  has  been  made 
of  light  emission  as  recorded  by  photodetectors,  pressure  or  tempera¬ 
ture  rise,  behavior  of  the  propellant  sample  after  interruption  of  the 


Price,  E.  W.,  and  others. "Theory  of  Ignition  of  Solid  Propel lan ta,M 
AMER  INST  AERON  ASTRONAUT  J,  Vol.  4,  No.  /  (July  1966),  pp.  1153-81. 
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ignition  stimulus,  and  physical  examination  of  the  specimen  after  quench- 
inf.  In  theoretical  models  the  attainment  of  a  given  surface  twrperature 
(often  called  the  Ignition  temperature) ,  rate  of  temperature  rise  or  of 
chealcal  heating,  or  mathematical  solution  of  the  equations  after  interrup¬ 
tion  of  external  heating  have  been  employed*  Only  United  work  has  been 
done  In  a  effort  to  correlate  the  different  Ignition  criteria,  but  it  Is 
known  that  the  definition  of  Ignition  is  important 


1-2  PURPOSE  OF  CURRENT  STUDT:  UNIFIED  IGNITION  THEORY 

While  there  has  been  reasonable  agreement  between  experimental  ig¬ 
nition  results  and  a  particular  theory  (probably  b>  purposeful  design). 

Lack  of  cocoon  and  compatible  assumptions  has  hindered  c  cm  par  Isons  among 
theories  and  between  theories  and  experiments.^  Hence,  the  same  experi¬ 
mental  data  have  seldom  if  ever  been  subjected  to  several  alternative 
theoretical  Interpretations.  One  purpose  of  this  Investigation  is  to 
facilitate  comparison  of  theories  and  experiment  by  constructing  a  model 
Including  chemical  processes  at  any  or  all  of  the  usually  selected  sites, 
but  otherwise  containing  compatible  and  consistent  ar sumption*.  Results 
of  solving  the  model  equations  are  not  expected  to  suostentlate  any  par¬ 
ticular  ignition  theory  to  the  exclusion  of  all  others;  Instead,  the  results 
will  probably  reveal  combinations  of  parameters  for  which  each  reaction 
site  (solid  phase,  surface,  gas  phase)  la  dominant  in  Ignition.  Many  of 
the  physical  quantities  Involved  are  net  accurately  known;  he  ice,  quanti¬ 
tative  testing  ci  the  theory  is  not  possible  at  present.  Nevertheless, 
the  unified  Ignition  model  should  enable  calculations  of  Important  sensi¬ 
tivity  relationships  between  parameters  and  solutions  as  well  as  meaningful 
comparison  with  results  from  appropriately  designed  experiments. 

A  second  important  our pose  of  the  study  Is  the  determination  of  the 
Importance  of  some  of  the  common  assumptions  which  are  often  made  to 
simplify  the  ma themat leal  treatment  of  the  equations.  Examples  are  in¬ 
clusion  or  exclusion  of  surface  regression,  gas  phase  convection  with 
consequent  blocking  of  heat  conduction  to  the  surface,  and  equality  of 
the  thermal  and  mass  dif fuslvltles  ( the  assumption  of  unity  Lewis  number). 


2.  PHYSICAL  DESCRIPTORS  of  problem 


Consider  the  transient  one-dimensional  system  in  Figure  1.  Initially, 
a  condensed  phase,  containing  fuel  and  oxidizer  ingredients,  occupies  the 
region  x  <  0  while  a  gas  phase,  consisting  of  oxidizer  and  Inert  ingre¬ 
dients,  fills  the  region  r  •  0.  In  general,  different  uniform  initial 
temperatures  are  permitted  for  each  phasr  At  time  zero  aa  ignition 
stimulus  consisting  of  a  time-dependent  flux  is  dl.ec^ed  at  Che  interface 
from  a  remote  position  In  the  gas  phase.  Condi  ions  remote  from  the  inter¬ 
face  are  maintained  at  their  Initial  values.  The  ensuing  response  of  the 
system  may  in  general  Include: 
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X-  0 


(GASIFICATION*  CHEMISTRY) 


FICUU  1.  Scbmtlc  Representation  of  Unified  Ignition  Model  . 


1.  Bulk  cbm  leal  react  Iona  confined  to  one  or  both  hooogeneou* 
phases*  In  the  solid  phase ,  the  reaction  occurs  either  betveen  fuel  and 
oxidizer  or  say  consist  of  a  pyrolytic  decomposition  of  the  solid  phase 
to  provide  a  reactive  fuel  and  oxidizer  which  participate  further  in  sur¬ 
face  or  gaseous  reactions*  In  the  gas  phase,  reactions  may  take  plica 
betveen  pyrolyzed  fuel  and  either  pyrolyzed  oxidizer  or  oxidizer  origi¬ 
nally  present  as  a  gas  phase  ingredient* 

2 .  Interface  processes,  consisting  cf  chemical  reactions  or  pyrolysis, 
in  which  two  phases  are  involved. 

3.  Motion  of  the  phase  boundary  owing  to  mass  transport  across  the 
bouadary.  In  this  model,  the  motion  Is  lndv  cd  solely  by  surface  processes 
either  pyrolytic  or  reactive  in  nature. 

4.  Gaseous  diffusion  of  all  species,  either  reactants  or  products. 
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5.  Absorption  of  th«  stimulation  radiant  flux  by  tbs  solid  phase  or 
by  evolved  gas**. 

6.  Hast  transfer  by  conduct loo  and  convection. 

In  tbs  following  mathematical  section,  the  appropriate  governing 
equations  are  developed  whose  solutions  describe  the  relationships  of  de¬ 
pendent  variables  (concentration  and  enthalpy  or  temperature)  to  indepen¬ 
dent  variables  (time  and  position). 


3.  KATKDUTICAL  FOWiULATlOH 


the  equations  which  describe  the  behavior  of  the  present  system  consist 
of  conservation  equations  and  constitutive  equations  combined  to  give  the 
governing  equations.  The  former  ere  general  and  Independent  of  the  nature 
of  materials  Involved;  the  latter  represent  empirically  observed  relation¬ 
ships  among  the  variables  appearing  in  the  conservation  equations  and  are 
required,  with  the  Initial  and  boundary  conditions,  to  provide  e  complete 
set  of  equations,  the  conservation  equations  may  be  further  classified  as 
field  or  Interface  equations. 

In  the  derivations  which  follow,  a  stationary  coordinate  system  relative 
to  the  observer  Is  established  at  the  phase  Interface.  Physically,  the  model 
behaves  as  If  material  In  the  solid  phase  approaches  the  boundary  from  the 
left  (see  Figure  1)  at  a  velocity  which  varies  with  time,*  crosses  the  phase 
boundary,  while  undergoing  substantial  physical  and  chemical  changes,  and 
flows  away  toward  the  right  with  a  velocity  which  Is  a  different  function  of 
time  and  position.  All  properties  are  assumed  constant  along  planes  parallel 
to  the  Interface,  so  that  the  model  Is  adequately  described  by  transient 
equations  of  one  space  dimension. 

3.1  COJCSBRVATIOH  EQUATIONS 

The  relationships  expressing  conservation  of  mass,  momentum  and  energy 
are  derived  snd  presented  in  numerous  standard  works  (for  example,  see  Foot¬ 
note  6).  In  the  current  study,  the  effects  of  gravity,  differential  exter¬ 
nal  force  fields  such  as  electrostatic  force  fields,  and  viscosity  are 
neglected;  snd  gas  velocities  are  assumed  to  be  low  enough  so  that  pressure 

*The  solid  phase  bulk  density  may  vary  due  to  differences  In  density 
among  the  reactants  and  solid  products,  but,  more  importantly,  due  to  sub¬ 
surface  evolution  of  gaseous  products.  Such  phenomena  are  neg  „•  ted  in  the 
present  study;  evolution  of  gas  is  presumed  to  occur  only  at  the  surface, 
and  solid  phase  reactions  are  assumed  io  occur  without  density  changes. 

For  these  reasons,  the  solid  phase  velocity  does  not  vary  with  position. 

6Bird,  R.  B.,  W.  E.  Stewart  and  E.  N.  Lightfoot.  Transport  Phenomena. 

Kev  York,  Wiley,  I960.  790  pp. 
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gradients  and  kinetic  energy  contributions  naad  not  b«  considered.  The 
usual  Euler  or  Kavler-Stc.\es  equation  of  notion  is  simplified  to  the  condi¬ 
tion  of  constant  pressure  (spacevlse).  The  nathenatical  description  of  the 
problem  then  Is  reduced  to  equations  representing  mass  and  energy  conserva¬ 
tion  for  each  bulk  phase  (field  equations)  and  the  Interface. 


3.1.1  Field  Equations 


3. 1.1.1 
stationary  elan 


Species  Conservation.  The  ns  as  balance  of  species  i  for  a 
rat  of  voluae  Is  given  by 


»(Prj)  »(pri°l) 


djr 


♦  at 


where  o  Is  the  density  of  the  medium,  end  Yj ,  Jtj  end  uj  ere  sees  free t ion, 
•rate  of  production,  end  velocity  of  species  i  relative  to  the  stationary 
coordinate  system.  By  summing  Eq.  1  over  mil  i  end  noting  that  ly*  -  1 
end  Skj  »  0 0  one  obtains  the  equation  for  overall  mesa  conservation 


3p  3(0v> 

St - 5 J— 

Equation  2  defines 


V  =  Ifj0j 


(2) 


(3) 


i 

i 


t 

t 

i 

! 

I 

! 

i 


as  the  mass  average  velocity,  usually  identified  aa  the  hydrodynamic 
(Pitot  tube)  velocity.  This  bar ycen trie  description  la  adopted  because 
of  its  formal  advantage  over  molar  or  other  descriptions  in  handling 
problems  dealing  with  simultaneous  convection  and  diffusion.  Ve  may  now 
write  Eq.  1  in  sn  alternate  fens  which  la  more  convenient  for  calculation 
purposes.  Ve  first  define 


r.  ?  u.  - 


(4) 


as  the  species  mass  velocity  relative  to  the  local  mass  average  velocity* 
The  product 


pVi 


=  h 


(S) 


represents  the  nass  diffuslonal  flux  of  species  i  relative  to  the  local 
average  velocity.  By  substituting  Eq.  4  and  5  Into  1  and  naming 
/  Eq.  2,  one  may  derive 


dr  3ri 

p  5T  H  p  5T  ♦  37- 


*r  +  ki 


(6) 
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as  the  expression  for  species  conservation. 

3. 1.1. 2  Energy  Conservation.  The  fundamental  derivation  of  energy 
conservation  is  in  terms  of  specific  internal  energy  of  a  stationary  volume 
element  and  leads  to 

*  tor)  -&-P&  <7, 


where  the  specific  internal  energy  V  is  made  op  of  thermal  and  chemical 
components. 

The  reduction  of  the  women tim  equation  to  p  ■  constant  makes  it  desir¬ 
able  to  eliminate  the  final  term,  containing  dv/lx,  from  Eq.  7.  This  may 
be  accomplished  by  making  use  of  Eq.  2  and  the  definition  of  enthalpy  (in¬ 
cluding  thermal  and  chemical). 


H  -  U  ♦  p/p 


The  final  result  is 


P 


da  . 
dt  z 


(8) 


(9) 


The  quantity  q  represents  heat  transferred  relative  to  the  mass  average 
velocic  /  and  in  general  Includes  heat  conduction  owing  to  temperature  and 
concentration  gradients,  transport  of  the  intrinsic  energy  of  diffusing 
species,  and  radiation.  In  principle,  photochemical  heat  release  could 
also  be  included  if  details  of  the  process  were  available  for  the  purpose 
of  establishing  the  effect  on  the  mass  conservation  equations.  The  final 
term,  dp/)t  is  what  remains  of  the  more  general  term 

5?  "  If  *  v  If  (10> 

vhrn  spacevise  pressure  variations  are  neglected.  In  the  following 
settlor**,  ;p/3t  is  also  ignored;  for  purposes  of  studying  dynamic  pressure 
effects  it  could  readily  be  reintroduced. 

3.1.2  Interface  Equations 


The  boundary  between  two  phases  is  physically  a  finite  region, 
through  which  the  properties  are  continuous  and  governed  by  the  same  con¬ 
servation  lavs  that  are  valid  fer  a  homogeneous  phase.  The  principal 
distinguishing  features  of  the  lnterfaclal  region  are  the  relatively  large 
spatial  derivatives  of  the  properties,  the  small  dimensions  compared  to 
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other  physical  dimensions  of  the  system,  and  s  Urge  change  in  relative 
importance  of  various  microscopic  force  terms  as  the  interfacial  tone  is 
traversed.  Examples  are  crystal  lattice  forces  and  binding  forces  of 
activated  chemical  coatplexes  associated  with  catalytic  reactions.  It  is 
convenient  In  the  current  discussion  to  bypass  the  details  of  such  phenom¬ 
ena  and  to  regard  the  interface  as  a  mathematical  surface  of  discontinuity 
joining  two  continuous  phases.  In  this  spirit,  it  Is  assumed  that  no 
accumulation  of  mass  or  energy  occurs  at  the  interface. 

3.1.2.1  Species  Conservation,  If  we  ignore  molecular  diffusion 
la  the  solid  phase,  the  mass  balance  of  species  i  at  the  surface  is 

mcim  *  *is  *  mgiM  #  *gin  Mil 


The  m  terms  represent  convective  mass  transport  owing  to  bulk  velocity  of 
the  phase  tdilch  arises  from  surface  chemistry  or  pyrolysis,  k  is  the  rate 
of  formation  of  species  i  in  surface  processes,  while  j  Is  the  gaseous  mass 
diffusion  rate.  Total  mass  conservation  is  obtained  by  summing  Eq.  11  over 
all  species : 


c  c 


-  p  V 

9  9 


(12) 


3. 1.2.2  Energy  Conservation.  The  Interface  equation  of  energy  con¬ 
servation  Is  readily  obtained  from  Eq.  9  by  dropping  the  ?p/3 1  term,  inte¬ 
grating  across  the  boundary,  and  setting  the  net  accumulation  to  eero.  The 
resulting  equation  la 


O  v  ff  ♦  <7  m  0  V  M  +  <7 

*C  c  c  ^ c  9  9  9  H9 


(U) 


which  indicates  that  convective  enthalpy  plus  heat  transport  relative  to 
hydrodynamic  velocity  Is  conserved  across  the  interface. 


3.2  CONSTITUTIVE  EQUATIONS 

Equations  6,  9,  11,  and  13  cannot  be  solved  as  they  stand  because, 
when  applied  over  all  indicated  species  and  relevant  phases  with  appropriate 
Initial  and  asymptotic  conditions,  there  are  more  unknown  quantities  than 
equations.  The  required  additional  equations  are  provided  by  the  transport 
relationships,  the  equation  of  state  (t^ken  for  convenience  as  the  perfect 
gas  lav),  kinetic  and  stoichiometric  expressions,  and  definitions. 
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3.2.1  Transport  Kelatlonshlps 

TIm  transport  relationship*  express  fluxea  In  terms  of  forcea  (grad¬ 
ient  a  of  physical  quantities  such  as  concentration,  enthalpy,  pressure, 
and  velocity).  A  fundamental  postulate  of  the  thermodynamics  of  irrever- 
aible  processes  is  that,  for  systems  not  far  moved  from  equilibrium,  the 
fluxes  are  linear  homogeneous  functions  of  the  gradients.  Expressions  of 
higher  degree  would  in  general  be  required  for  large  departures  of  the 
system  from  equilibrium;  however,  little  work  has  been  done  on  the  higher 
order  theories*  Experimental  results  have  confirmed  the  adequacy  of  the 
linear  relation ships  over  a  vide  range  of  conditions* 

In  the  absence  of  viscous  forces,  the  fluxes  to  be  considered  ere 
mats  and  energy.  The  forces  Involved  are  related  to  gradients  of  pressure, 
temperature,  a^d  concentration.  In  the  present  analysis,  only  the  direct 
effects  are  considered;  the  effects  of  pressure  gradient  on  mass  flux,  of 
temperature  gradient  on  mass  flux  (thermal  diffusion  or  Soret  effect),  of 
external  forces  on  mass  flux  (forced  diffusion),  and  of  concentration 
gradient  on  energy  flux  (Du four  effect)  are  neglected.  The  Soret  effect 
is  the  most  significant  under  ignition  conditions  but,  even  so,  it 
probably  contributes  a  maximum  of  only  10X  of  the  total  mass  flux.  Vlth 
these  simplifications,  the  energy  and  mass  fluxes  for  a  multi-component 
system  are  given  by 


(14) 


and 


n 


(IS) 


In  Eq.  14,  the  three  terms  on  *he  right  express  energy  flux  due  to  normal 
thermal  conduction,  diffusing  species,  and  ndf  itlon,  all  relative  to  the 
barycentrlc  velocity  v.  Equation  IS  is  obtained7  from  a  more  general  ex* 
press  ion  by  assuming  the  perfect  gas  lav.  It  is  seen  that  for  a  multi¬ 
component  system,  the  mass  diffusion  fluxes  of  the  individual  species  are 
complicated  functions  of  the  concentration  gradients  of  all  the  species. 
The  multicomponent  diffusion  coefficients,  Djj,  are  generally  functions 
of  all  the  concentrations.  Equation  IS  may  be  expressed  in  different 
form^by: 


me  TT  5618,  T* rt  1 


n 


(U) 


Equation  16  represents  a  sat  of  simultaneous  equation*  relating  the  mat* 
diffusion  fluxes*  the  mat*  concentration*  and  the  mat*  concentration  gradient* 
7r.  An  advantage  of  E q.  16  1*  that  the  binary  diffusion  coefficients  V  are 
nearly  Insensitive  to  concentration;  disadvantages  are  that  the  fluxes  are 
expressed  implicitly  and  that  n(n-l)/2  binary  coefficients  must  be  specified 
to  solve  the  probla.  This  latter  requirement  greatly  Increases  the  number 
of  quantities  required  for  a  parameter  study;  other  Investigations8  have 
shown  that  often  the  principal  trends  of  results  may  be  adequately  calculated 
by  assuming  all  binary  coefficients  to  be  equal.  With  this  asswptlon  Eq*  15 
or  16  may  be  reduced  to  the  form  used  in  this  report 


(17) 


triilch  is  knovn  as  Flck’s  first  lav  of  diffusion.  Ve  note  that  the  deriva¬ 
tion  of  Eq.  17  does  not  depend  upon  the  frequently  made  assumption  of 
equal  molecular  weights  of  all  species. 

3.2.2  Kinetics  and  Stoichiometry  . 

The  conservation  equations  contain  terns  representing  rates  of  trans*> 
formation  of  chemical  species  (*j,  Jtj,  for  bulk  and  surface  reactions)  and 
rates  of  mass  evolution  (manifested  by  the  convective  terns  in  the  interface 
equations).  Chemical  transformations  are  complicated  by  multiple  path* 
opposing*  or  consecutive  reactions*  details  of  vh  ;h  are  not  fully  knovn  In 
complex  combustion  processes.  The  assumptions  regarding  the  kinetic  lavs 
idilch  must  be  made  In  order  to  achieve  tractablllty  probably  reflect  the 
greatest  art~  of  uncertainty  in  the  formulation  of  any  Ignition  or  combustion 
theory.  It  la  ;onvenient  to  assume  that  the  rate  controlling  step  at  a  given 
site  In  a  kinetic  process  consists  of  a  simple*  single  step.  Irreversible 
reaction  with  rate  folloving  an  Arrhenius  dependence  on  local  temperature. 

The  stoichiometry  of  the  overall  chemical  change  vlll  not  correspond  to  that 
given  by  the  reaction  mechanism  If  other  faster  reactions  occur  Involving 
the#  same  components.  We  nov  consider  the  kinetic  expressions  for  the  changes 
appearing  In  the  conservation  equations.  The  processes  to  be  considered  were 
enumerated  In  Section  2. 

3.2.2. 1  Gas  Phase.  The  two  gas  phase  reactions  Included  in  the  formu¬ 
lation  are; 


Kendall  R.  M. ,  R.  A.,  Rlndai  and  E.  P.  Bartlett.  MA  Multicomponent 
Boundary  Layer  Chemically  Coupled  to  an  Ablating  Surface,”  AKER  INST  AERON 
ASTRONAUT  J,  Vol.  5,  No.  6  (June  1967),  pp.  1063-71. 
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* 

» 


Reaction 

7t 

iri  *  b2toi  -  b7i?i 

(18) 

Reaction 

$t 

iri  *  b4to i  -  b9in 

(29) 

Equation  18  represents  react ion  of  gas  phase  fuel  (pyrolyzed  from  nolid) 
with  oxidizer  pyrolyzed  fro*  the  solid  phare.  The  h's  ere  molar  stoichio¬ 
metric  coefficients  relative  to  luel  taken  aft  one.  Mechanistic  reaction 
orders  with  respect  to  fuel  and  oxidizer  are  taken  as  F7  and  07 9  respectively. 
Equation  19*  employing  analogous  notation*  pertains  to  the  reaction  o*  pyro- 
lyced  fuel  with  oxidizer  originally  present  in  the  gas  phase;  reaction  orders 
are  denoted  as  F8  and  08.  Mass  balance  considerations  lead  to: 


UJ  *  V?2  m  u7b7 

(30) 

V1  *  m 

(21) 

wb.r.  th.  u's  draott  aoltcultr  vtlghu,  v<  referring  co  th.  molyttd  kuwui 
fool.  Applying  th.  lw  of  mu.  .ctlon  of  chemical  kinetics*  to  Reactions  7 
and  8,  w*  nay  represent  th.  um  .cure,  terms  (cf.,  Eq.  6)  as: 

Gas  phase  fuelt  k ^  -  *wj(r7  *  rg) 

(22) 

Gas  phase  oxidizer 

(pyrolyzed  from  solid)  t  k ^ 

(23) 

Gas  phase  oxidizer 

(initially  in  gas)  t  k*  »  -u^b^r^ 

(24) 

Product  of  reaction  7:  k?  - 

(25) 

Product  of  reaction  f: 

(26) 

where 

r7  *  V7 iVl'Vl)r7  (W/7  «P('V*T) 

(27) 

r8  -  ivs^)08  «p(-V*^ 

(28) 

®  VUllans,  F.  A.  Cofrbuaticn  Theory.  Palo  Alto,  Ca 
Pub.  Co.,  1985,  AA7  pp. 

.,  Addlson-Weslcy 
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and  where  v  and  F  represert  pre-exponential  factor  and  activation  energy 
of  the  appropriate  reaction.  Negative  signs  in  the  expressions  for  k  In¬ 
dicate  disappearance  of  the  relevant  species.  Species  3,  5,  6,  and  9  are 
not  generated  by  gaseous  rsictlon;  consequently,  their  corresponding  Jt*s 
are  aero. 


3.2. 2.2  Solid  Phase.  Reaction  kinetics  In  the  solid  phase  are 
no  re  complicated  than  In  the  gas  phase.  The  components  of  a  solid  are 
not  completely  mobile  as  Is  the  case  for  a  gas;  Instead,  Individual 
molecules  are  constrained  to  move  relatively  small  distance  within  crystal 
lattices  or  across  interfaces  separating  different  constituents.  Since 
each  molecule  does  not  ha^  s  ejual  probability  of  collision  with  every  other 
molecule,  the  usual  method  used  in  deriving  rate  expressions  for  gases 
cannot  be  employed.  Therefore,  solid  phase  reactions  may  be  expected  to 
depend  upon  first  order  decompositions  or  upon  reactions  which  occur  at 
internal  interfaces.  In  the  latter  case,  the  question  of  reaction  rate 
cannot  be  resolved  without  reference  to  geometrical  considerations.  Even 
the  extent  to  which  the  reaction  can  proceed  is  limited  by  the  diffusion 
rates  of  the  solid  reactants  through  films  of  products  built  up  at  the 
Interfaces.  It  Is  clear  from  discussions  nresented  lathe  literature10*11 
that  no  universally  applicable  rate  theory  Is  avsllable  and  that  topocheni- 
cal  considerations  would  be  required  to  provide  an  accurate  description  of 
the  solid  phase  process.  However,  the  genersl  aspects  of  solid  phase  re¬ 
actions  can  be  expressed  by  assuming  reaction  orders  ranging  from  aero  to 
two.  Moreover,  it  has  been  shown1*  that.  Insofar  as  Ignition  Is  concerned, 
the  thermal  history  of  a  solid  undergoing  bulk  reaction  is  relatively  In¬ 
sensitive  to  the  kinetic  order  assumed.  For  simplicity  In  the  present 
analysis,  the  following  assumptions  are  made:  (1)  the  kinetics  Is  sero 
order,  and  (2)  the  solid  phase  composition  can  be  characterised  by  a 
fuel-to-oxidlxer  ratio  (even  for  homogeneous  double-base  propellants). 
Assumption  (1)  requires  that  careful  account  be  kapt  of  the  amount  of  re¬ 
actant  consumed  to  prevent  the  physically  unrealistic  and  numerically 
erroneous  Introduction  of  negative  concentrations  during  computations. 

The  single  solid  phase  reaction  is  assumed  to  be 

Reaction  10 1  IF]  ♦  *nlOl  -  (**) 

with  the  product  [F]  considered  ss  inert,  or  the  veverse  reaction: 

b12lpI  ■  trj  *  bntoi  (jo) 


Corde* ,  Herman  F.  "The  Preexponent 1*1  Factor*  for  Solid-State  Thermal 
Decoapositlon,"  J  CHW  PHYS,  Vol.  72,  No.  6  (June  1968),  pp.  2185-89. 

H  Caraer.  W.  E.,  ed.  Chorietry  of  the  Solid  State.  London,  Buttervorths, 
1955. 

12  - 

Bradley,  H.  H..  Jr.  HTheory  of  Ignition  of  a  Reactive  Solid  by 
Com t ant  Energy  Flux,"  COMBOS  SCI  AJ©  TECH,  Vol.  2,  No.  1  (August  1970), 

pp.  11-20. 
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representing  the  bulk  phase  dissociation  Into  active  fuel  and  oxidizer. 
Hass  balance  yields 

U10  *  UJJbJJ  "  UJ 2*12  (U> 

and  the  source  teres  for  reaction  10  (Eq.  6)  are  given  by 


kio  •  ’uiorio 


kJJ  *  ~vuburio 


k12  m  U12»12r 


10 


(32) 

(33) 

(34) 


The  rat*  rj0  la  given  by 


rio  *  vj 


OS) 


where  Is  the  activation  energy  of  the  solid  reaction.  The  reverse 
process  may  be  readily  simulated  by  changing  the  signs  of  the  expressions 
for  the  k'*  in  Eq.  32-34  or  of  vj  in  Bq.  35. 

3.2.2. 3  Interfacial  Processes.  The  boundary  between  solid  and 
gas  phases  may  be  the  site  of  chemical  transformations  or  of  phase  changes 
Involving  either  simple  gasification  or  dissociation.  Rates  of  chemical 
transformation  at  the  Interface  depend  upon  adsorption  and  desorption  pro¬ 
cesses.  Details  of  surface  phenomena  are  discussed  In  the  literature? 

(see  p.  382  et.  seq.);  the  results  may  be  summarized  by  the  statement  that 
the  apparent  order  of  the  reaction  can  vary  between  zero  and  the  true 
kinetic  order  and  nay  be  of  order  minus  one  relative  to  the  reaction  pro¬ 
duct.  The  latter  occurs  when  the  product  is  strongly  absorbed,  thereby 
Inhibiting  the  adsorption  of  additional  reactants.  The  surface  reactions 
considered  in  this  analysis  are:  (1)  reaction  of  pyrolyzed  oxidizer  with 
solid  fuel  and  (2)  reaction  of  original  oxidizer  from  the  gas  phase  with 
solid  fuel.  They  are  represented,  respectively,  by  the  following  stoichio¬ 
metries: 


Reaction  5; 

ir) 

♦  b'2(0) 

-  b5(P 1 

(36) 

Reaction  6: 

tr) 

♦  b'4ioi 

-  b6(P) 

(37) 

Corresponding  to  reactions  5  and  6,  we  have 
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PJ0  ♦  '  y5b5 

U10  *  Vi  "  V> 6 


(38) 


(39) 


The  fuel  mees  consumption  races  for  Che  above  are  given  by: 


*5  '  pcr!0V5  (P9r/M^55  "P'-V". i 


(40) 


iw-r  -p{-v".) 


(4iJ 


It  is  assumed  that  the  surface  rate*  are  first  order  with  respect  to  the 
fuel  (proportional  to  fuel  surface  area)  and  of  order  S5  and  S6  with  re¬ 
spect  to  the  pyrolyzed  or  initial  oxidizer. 

Ingredients  subject  to  surface  pyrolysis  Include  initial  solid  fuel 
and  oxidizer  together  with  their  solid  phase  reaction  products  *  if  any. 
The  rates  of  pyrolysis  may  be  controlled  by  an  A  rhenlus  temperature 
expression  or  may  proceed  by  an  equilibrium  vaporization  process  in  which 
the  gas  phase  equillbrlua  vapor  pressure,  hence  sole  or  mass  fraction,  is 
related  to  total  applied  pressure  and  surface  temperature.  Combinations 
of  the  two  processes  lead  to  intermediate  results.  The  assumption  of 
surface  rate  controlling  processes  leads  to  more  loosely  coupled  boundary 
conditions  and  is  adopted  here  for  reasons  of  mathematical  simplicity. 

A  brief  description  for  the  case  where  at  least  one  equillbrlua  process 
is  involved  is  presented  in  Appendix  A.  The  mass  fluxes  of  the  solid 
phase  fuel,  oxidizer,  and  product  are  given  for  i  -  10,  11,  and  12,  re¬ 
spectively,  by: 


P(- ) 


(42) 


ai 


where  fractional  surface  area  of  ingredient  i  has  been  set  equal  to  its 
mass  fraction.  The  values  of  YC£S  depend  in  general  upon  the  initial 
composition  of  the  solid  phase  and  the  extent  of  solid  phase  reaction  and 
so  depend  upon  a  solution  of  the  solid  phase  equations. 

The  rigorous  application  of  Eq.  39-42  involves  consideration  of  the 
changing  topology  of  a  heterogeneous  surface.  The  complication  is  brought 
about  by  the  fact  that,  unless  the  kinetic  parameters  and  surface  areas 
of  the  different  components  are  uniquely  related,  each  component  oust  burn 
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at  a  different  linear  rate,  a  conclusion  which  violates  the  assiaptlon  of 
one  dimensionality*  In  steady  state  burning.  It  has  been  speculated  that 
•ach  component  burns  at  a  different  temperature^  bur  this  explanation  Is 
not  universally  accepted*  Actually,  it  la  known  that  In  steady  state,  each 
component  burns  at  a  different  rate,  with  the  surface  topology  adjusting 
Itself  to  give  areas  necessary  to  accommodate  equilibrium  mass  ratios 
typical  of  the  solid*  Further  discussion  to  resolve  the  problem  does  not 
seem  appropriate  at  present*  The  simplifying  assumption  Is  therefore  made 
that  each  component  of  the  propellant  regresses  at  a  rate  consistent  with 
the  kinetic  equation;  the  degree  of  departure  of  the  surface  from  planar 
can  then  be  partially  assessed  by  examining  the  total  regression  of  each 
component  at  the  time  of  ignition  and  comparing  differences  with  reasonable 
particle  sixes.  If  such  a  test  reveals  that  the  surface  has  indeed  become 
two-dimensional,  then  a  two-dimensional  analysis  must  be  made  or  the  alter - 
mate  assumption  adopted  that  the  regression  rate  of  the  fastest  burning 
component  governs  the  overall  regression  rate.  In  addition.  It  Is  asswed 
for  convenience  that  no  solid,  unreacted  material  Is  ejected  from  the  sur¬ 
face  to  contribute  a  thermal  effect  In  the  gas  phase*  It  should  be  noted 
that  in  no  one-dimensional  analysis  Is  It  possible  to  consider  mechanisti¬ 
cally  the  changing  surface  topology  together  with  associated  changes  In 
surface  areas  of  the  components* 

3.2.3  Radiant  Energy  Transfer 

Except  under  hypergollc  conditions  at  ambient  temperatures,  some 
sort  of  external  snergy  source  serves  as  the  ignition  stimulus  In  Ignition 
experiments*  Ordinarily  this  consists  of  a  radiation  source  (laser,  xenon 
arc,  or  carbon  arc)  focused  outo  the  surfree  of  the  sample  to  be  tested* 

The  fractions  of  this  flux  which  are  absorbed  and  reflected  depend  upon 
the  optical  properties  of  the  surface  and  of  the  bulk  solid  phase.  Absorp¬ 
tion  nay  be  confined  to  the  surface  for  a  completely  opa'jue  material  or 
one  with  an  absorbing  agent  applied  to  the  surface*  On  the  other  hand, 
the  energy  may  be  deposited  Instantaneously  In  depth  for  a  t  r  anal  use  ent 
naterlal,  with  attenuation  being  adequately  described  by  an  exponential 
relationship  between  flux  density  and  distance  (Beer9*  Lav).  With  the 
passage  of  tine,  gases  evolved  from  the  surface  processes  commence  to 
attenuate  the  incident  energy  bean.  Again,  If  re-radlatlon  by  the  gaaaa 
la  not  appreciable,  exponential  decay  of  intensity  is  a  reasonable  approx¬ 
imation.  This  attenuation  of  the  Incoming  st 'cuius  by  the  evolved  gases 
adds  energy  to  the  gas  phase  while  denying  it  to  the  heated  surface*  As 
the  temperature  of  the  gas  phase  rises,  re-radlatlon  by  the  gases  also 
begins  to  Intervene.  All  these  phenomena  are  complicated  by  being  wavelength 
and  species  dependent,  vita  gaseous  absorption  of  reflected  energy  also 
playing  a  role  of  some  importance.  The  principal  affect  Investigated  here 

**  Schultz,  R.,  snd  L.  Green,  Jr.  Third  AGAFD  Colloquium,  Pergamon 
Press,  l9S8.  "Studies  of  the  Decomposition  Mechanism,  Erosive  Burning, 

So nance,  and  Resonance,  for  Solid  Composite  Propellants,  In  Combustion 
and  Propulsion,*9  pp*  367-420. 
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Is  the  absorption  of  energy  by  the  solid.  It  la  assumed  that  Incident  flux 
density  is  spatially  and  temporally  constant,  and  that  the  surface  reflect¬ 
ivity  remains  constant  during  the  heating  phase.  A  simplified  analysis  of 
the  inclusion  of  gas  phase  energy  attenuation  is  presented  in  Appendix  B. 

3.2.4  Relationship  of  Enthalpy  to  Temperature 

In  the  foregoing  sections,  the  energy  equation  has  been  presented 
In  tens  of  total  enthalpy  while  reaction  rates  involve  temperature,  which 
is  related  to  the  thermal  enthalpy.  The  general  relationship  between  the 
temperature  and  total  enthalpy  is 


(43) 


where  Af  is  the  chemical  enthalpy  of  species  i  at  some  arbitrary  refer¬ 
ence  teaperature  Tr.  One  convention,  adopted  here,  is  to  assign  zero  en¬ 
thalpy  to  the  elements  at  Tr,  so  that  hj  is  the  heat  attending  the 
formation  at  Tr  of  species  i  from  the  elements. 

3.3  GOVERNING  EQUATIONS 

The  relationships  presented  in  Section  3.1  and  3.2  describe  the 
phenomenon  of  solid  propellant  ignition  completely,  subject  to  the 
llmltat/^ns  mentioned.  Owing  to  the  nonlinear  nature  of  the  equations, 
exact  analytic  solutions  have  net  been  found;  recourse  eust  be  to  numeri¬ 
cal  methods  or  to  one  of  several  available  approximate  analytic  techniques. 

A  desirable  resolution  of  the  probleza  would  be  the  application  of  both 
approaches  and  comparison  of  results;  however,  only  the  numerical  method 
is  to  be  employed  In  the  current  report. 

While  numerical  analogs  may  be  developed  for  the  equations  of 
Sections  3.1  and  3.2  as  they  stand,  it  Is  more  convenient  to  combine  some 
of  them,  to  apply  certain  mathematical  transformations,  and  to  make  ad¬ 
ditional  simplifications  In  the  Interest  of  reducing  the  number  of  physi¬ 
cal  quan^lt'es  to  be  subjected  to  a  parameter  study. 

3.3.1  Species  Conservation 

Equations  6  and  17  may  be  combined  directly  to  give  a  single  equation 
describing  the  ti&e  and  space  distribution  of  each  species.  For  the  gas 
phase  (i  -l,  2  ...  9;  x  >_0)  the  result  is 


(44) 
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In  tha  solid  phasa  (i  -  10,  11,  12  x  <  0)  where  the  aau  diffusion  Is 
negligible,  tbs  species  halar^es  ars 


P  5 

vc  dt 


-  kt 


(4S) 


Equation  11,  applied  to  each  epee  lea  ,  yields  tlie  following  Interface  con¬ 
servation  lavs: 


3f. 


**U  m10  •  Vi  '  Pcp  7T 


146J 


5r-» 

Solid  oxidizer:  •  t  f,  -  pt-r"*  : - 

JI  *  2  g  Tx  MJ0  2"5 

3r 

Solid  product:  »„  *  »  T,  -  P  ®  *5 — 

12  a  3  g  lx 


(47/ 


i 48 ; 


3y.  w. 

Caseous  oxidizer:  0  •  »  r  -  p  P  -r —  ♦  - t's. 

a  4  vg  lx  u10  4  6 


(49) 


Product  of  surface 

reaction  of  fuel  with  Bf-  U- 

py  roly zed  oxidizer :  0  -  m  Y-  -  p  V  - - b  (50J 

s  5  Kg  dx  Uf/1  >5 


Product  of  surface 

reaction  of  foei  with 
gaseous  oxidizer:  0  m 


(SI) 


Product  of  gaseous 

reaction  of  fuel  with 
pyrolyzed  oxidizer:  0  «  ms*j 


Product  of  gaseous 

reaction  of  fuel  with 
gaseous  oxidizer:  0  m  ms*$  " 

3rg 

Inert  gas:  0  ■  ■  I  *  p  P  -t? — 
*  s  9  g  dx 


3r7 

-  V  IT  (S2) 

3yfl 

o  o  _ - 

Pg  3*  rsi; 

(S4) 
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Equation  (46)  deserves  special  cotst.  Although  the  total  mss  flux  of 
fuel  contributing  to  s^is  in  Eq.  11  includes  *5,  n$,  and  m^,  the  reaction 
terns  1*5  and  are  exactly  cancelled  by  Jt^s.  The  quantity  n*  In  Eq.  46*54 
represent  a  the  sua  of  all  case  transferred  across  the  boundary  so  that 
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11 


12 


(SS) 


The  gaseous  co erection  ten  in  each  surface  conservation  equation  is  rep* 
resented  ay  %rj,  where  the  on^d leans ional  aasuaptlon  of  complete  mixing 
of  gaseous  species  at  the  surface  has  been  made. 

3.3.2  Enerer  conserratlon 

If  we  assime  specific  heats  are  independent  of  temperature,  it  is 
possible  to  combine  Eq.  6,  9,  14,  17,  and  43  to  obtain  the  following  equa* 
tlon  for  thermal  enthalpy  h; 


I 

XX 

XXX  17 

&|c^[s7*  '“-"“j  ir] 

(56) 

where 

“l  *  Ci(r  -  Tr) 

(57) 

1 

> 

9 

(58) 

Im  -  PC  D/A 

av 

(59) 

C  *  Ec  f 
ar  ^ixi 

(60) 

The  terms  in  Eq.  56  have  the  following  physical  interpretations : 

Thermal  energy  accumulation  in  a  moving  volume  element 

II  Thermal  energy  added  to  volume  element  by  differential  processes 

III  Chemical  energy  generated,  alternately  represented  by  V^OjTj* 


where  Pj  and  rj  are  tht  exothermicity  and  rate  of  the  j th  reaction 
with  Oj  evaluated  at  the  reference  temperature.  The  factor  vr  la 
the  molecular  weight  of  the  ingredient  with  respect  to  which  heat 
effects  are  computed. 
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IV  Radiant  energy  deposited  in  depth  for  the  eolld  phase.  The  anal* 
ogous  ten  la  discussed  in  Appendix  i  for  the  gas  phase. 

The  portion  of  the  diffuslotul  tern  containing  Le-I  as  a  factor  accounts 
for  two  phenomena: 

(1)  difference  In  dlffusional  fates  of  mass  and  energy 
,  (2)  difference  in  specific  heats  of  reactants  and  products 

It  is  Important  to  note  that  In  the  Chemical  energy  term,  the  reaction  rate 
Jtj  la  taken  at  the  Instantaneous  local  temperature  while  the  heat  effect 
la  related  to  the  reference  temperature.  The  fact  that  heat  of  reaction 
la  generally  temperature  dependent  is  associated  with  the  tern  containing 
I*-2  as  a  factor.  In  the  Interest  of  mathematical  simplicity,  we  shall 
simplify  Eq.  56  by  the  single  assumption  of  equal  specific  heats  for  all 
species.  This  is  tantamount  to  assuming  a  heat  of  reaction  which  is  large 
compared  with  sensible  heats  over  the  range  of  temperature  involved.  The 
assumption  of  unity  Levis  number  is  then  unnecessary  as  the  expression 
^*•2  vanishes.  Moreover,  Cav  becomes  constant,  thereby  further 

simplifying  the  ensuing  numerical  procedures.  The  major  effect  of  Levis 
oomber  still  appears  in  Eq.  44  an d  in  the  mass  interface  conditions  vhen 
pJP  !•  replaced  by  its  equivalent  from  Eq.  59.  The  effects  of  different 


specific  heats  My  be  studied.  If  desired  by  using  Eq.  56  as  s  point  of 
departure  and  considering  as  a  function  of  x,  thereby  Introducing  local 
values  of  and  3r^/3x  into  Eq.  56. 

The  boundary  condition  for  the  energy  equation  la  formulated  by  combining 
Eq.  13,  14,  17,  and  43,  with  the  assumption  of  equal  specific  heats  for  both 
phases.  This  latter  step  is  equivalent  to  assuming  that  surface  heat  effects 
(chemical  reaction  and  phase  changes)  are  large  compared  vith  changes  of 
sensible  enthalpy.  The  final  equation  is 


(61) 


Although  apparently  inconsistent  vith  the  above  assumption,  unequal  specific 
heats  are  retained  for  the  terms  of  the  two  phases  in  Eq.  61  because  unequal 
specific  heats  represent  only  a  second  order  effect  on  qs,  the  surface  heat 
ralease  term,  whereas  the  specific  heat  is  a  first  order  effect  in  the  tents 
for  conductive  energy  transport.  The  term  qT  is  the  radiant  hast  absorbed  at 
the  mathematical  interface.  This  thermal  contribution  is  included  either 
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la  Eq.  56  vtwa  0  1*  flnit*  or  la  tq.  oi  vaaa  V  la  Infinite,  but  not  la 
both.*  Tbn  surface  has  tins  qa  la  glv«n  by 


'  m5°5  *  m6°6  *  Ml(PlO  *  mU°U  4  mlfl2 


(62) 


in  which  ths  firs  Q'm  arc  the  exothermic  it  las  accompanying  the  various 
surface  changes  (see  Section  3.2. 2*3) . 

3*3*3  Auxiliary  Conditions 

Tb  complete  the  ays  tee  of  field  and  interface  equations  *  initial  con* 
dittoes  and  conditions  in  each  phase  reeote  froe  the  lnterf see  are  pre¬ 
scribed. 

Initial  Conditions: 

Solid:  Gas: 

^  (x$o)  m  r^o  ^  (x*o)  m 

fuel:  Y10  (x,o)  -  Yfc  Oxidiser:  r4  (x,o)  -  1  -  Y^ 

Oxiditer:  Ylx  (x.o)  •  1  -  Yfc  Inert:  Yg  (x^>)  -  1  - 

Product:  Yi2  (x&)  m  0  All  other  YX(i  •  1-9)  *  0 

laeote  Conditions:  Gradient!  of  all  quantities  yanish  at  x  »  — •  for  the 
solid  phase  and  at  x  -  «•  for  the  gas  pha~*. 


*  Mathematically,  as  0  approaches  infinity,  the  behavior  of  the  solo* 
tlon  of  the  equations  is  as  if  the  radiation  tern  disappeared  from  the 
field  equation  and  appeared  in  the  boundary  condition.  Finite  difference 
aethods  fall  to  realize  this,  however,  and  the  dual  treatment  becomes 
necessary  because  of  limitations  of  computers  to  the  representation  of 
finite  quantities. 


\ 
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6.  MANIPULATION  OF  EQUATIONS 

The  foregoing  sections  have  been  devoted  to  the  development  of  the 
fundamental  equations  of  a  unified  Ignition  theory.  Presented  In  their 
physical  fora  they  are  valuable  In  understanding  the  details  of  the 
physical  processes  Involved;  however,  appropriate  modifications  are  needed 
to  convert  the  equations  to  a  more  convenient  form  for  computational  use. 

In  general,  these  include  non-dimensions 11 ration  through  the  use  of  charac¬ 
teristic  quantities  and  coordinate  transformations  such  as  the  use  of 
similarity  variables  and  stretching  parameters.  The  equations  at  hand 
reprssent  such  a  diversity  of  physical  phenomena  that  there  Is  no  single 
set  of  characteristic  quantities  which  apply  In  general.  The  resolution 
of  this  problem  will  be  to  delay  the  normalization  procedure  until  a  parti- 
cola*'  subcase  of  the  unified  theory  has  been  chosen  for  study.  In  this  way 
ve  retain  flexibility  while  avoiding  the  choice  of  a  characteristic  quantity 
whose  value  may  approach  a  limit  of  rero  or  Infinity.  Similarity  transfor¬ 
mations  have  not  yet  been  found  for  the  type  of  equation  Involving  an  ex¬ 
ponential  nonlinearity  In  the  dependent  variable.  This  section  will  present 
the  two  coordinate  transformations  used  to  convert  the  equations  for  compu¬ 
tational  use. 


4.1  HDVARTH  TRANSFORMATION 

The  substantial  or  total  derivatives  which  appear  In  Eq.  44,  45,  and  56 
contain  the  local  phase  velocity  Implicitly  (see  Eq.  6  and  9).  Owing  to  the 
constant  pressure  approximation  of  the  momentum  equation,  only  the  global 
continuity  equation  remains  as  an  independent  relation  containing  the  local 
velocity  Let  us  define 


(63) 


with  p0  taken  as  some  reference  (here  initial)  value  of  p  whence 


(64) 


and 


a 

dM  0 


(65) 


26 


WK  TT  561$,  Fart  1 


or 


} 

Ti 


p  a 


Then  oslng  tb.’ following  standard  rule*  of  partial  dlffaraotlation: 


(67) 


(*•> 


wa  say  darltra 


d_ 

dt 


■  £  (&),  *  (n). 


(«» 


Application  of  Eq.  (66)  and  (69)  to  Eq.  (46),  (4$)  and  ($6)  provldaa  tba 
following: 


,ri  .  ,rt .  j  j|„ •>„  ,*i 

po 

3ft  _  **  3ft  A  1  3  (p\  3ft  \  A  ur  6  *  I  8  A 

at  •  -  r  5?  ♦  77  5?  (t  s*r  T  *  T  qr  **P  (po  T  *) 


f7o; 


f7I/ 


Th«  latcrfici  Eq.  46-54  are  not  greatly  changed  except  that  the  diffusion 
terms  becc 


p  V arf/a* 


In  Eq.  71  the  term  involving  Lewis  number  has  been  eliminated  for  reasons 
previously  discussed  (see  Section  3.3.2) .  The  final  term  of  Eq.  71  vould 
appear  only  for  the  solid  phase  and  for  finite  6.  Equation  70  applies  to 
the  solid  phase  with  V  sst  to  zero. 

Results  of  using  the  Bovarth  transformation  era: 

(1)  The  global  equation  of  continuity  Is  satisfied  identically  in  the 
transform  space.  As  a  result  the  set  of  Eq.  70  contains  one  redundant 
equation  of  each  phase.  This  redundancy  is  retained  in  the  numerical  an¬ 
alysis  to  provide  an  accuracy  check. 
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(2)  Th.  coefficient  of  th«  convective  t«rm  (>/!♦)  Is  not  t  function 
of  position  but  depends  only  upon  tbo  lnotantanooua  velus  of  tba  total 
rata  of  naaa  tranafar  across  tba  lntcrfac*. 

2  (3)  Tba  cosffldants  appearing  In  tha  diffusion  tarns,  i.r.. ,  pX/c  and 

p  V ,  ara  nors  naarly  lndapandant  of  tMparatura  for  a  (Ivan  pr as sura  than 
ara  tbalr  counterparts  In  the  untransfornad  aquations;  hanca,  thay  nay  ba 
sored  outside  tba  partial  differential  operator. 

A  disadvantage  of  Bovarth  transformation  is  tba  need  to  par ft  m  tba 
Inverse  (Integral)  transfornatlon 


*  » 


fk  f*—!i _ 

a  J0 


(72) 


in  order  to  recover  the  true  distance  scale.  If  the  surface  conditions  are 
all  that  are  required,  this  extra  step  is  eliminated. 


4.2  TRANSFORMATION  FROM  IK?  I  KITE  TO  FINITE  DOMAIN 


The  extent  of  the  space  scale  ranges  frou  minus  to  plus  infinity.  Bov- 
ever,  computer  simulation  by  finite  differences  can  deal  only  with  finite 
quantities.  One  way  out  of  the  difficulty  is  to  repeat  the  cosqxstional 
effort  for  several  large  finite  values  for  the  distance  of  the  rsaote 
boundaries  until  the  solutions  converge. 1*  The  technique  teployed  here  is 
to  use  an  exponential  transformation  of  the  distance  scale  which  has  the 
added  feature  of  providing  closer  spacing  of  space  nodes  near  the  interface 
where  gradients  are  larger.  The  transformation  relationships  are 


♦c  ■  *2  ♦ exp  Me) 


(73) 


(74) 


Me) 


for  the  solid  phase  (-1  <  $  <  0) 

c 


(75) 


Kav.l  Weapons  Center.  Theory  of  Ignition  of  A  Reactive  Solid  by  Constant 
Energy  F^a:,  by  H.  H.  Bradley,  Jr.,  China  Lake,  Calif.,  NWC, November  1968. 

44  pp.  (NWC  TP  4618). 
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and 


v J  •  •«’  K*») 

JF  ‘  0»(J’*«)  5C 


(7«i 


f77; 


V  •  -  3^] 


(7i) 


(or  tba  |u  pbaaa  (0  <  6^  <  27. 

Application  of  Eq.  73-76  yields  tba  (Inal  transformed  equations  to 
ba  solved  numerically: 


J*i  ,  ***,  Hi 

IT  -  -TT  '  W  *  vi 


H 


with  Initial  conditions 


*i<*'  °>  -  *io 


M  -  1-14) 


(1  -  2-247 


and  remote  boundary  conditions 

0  1  *  l->r  17  pas 

'  i  •  10-12,  14  ao2id 

Interface  boundary  conditions  for  species  are  given  by 


(79) 


(80) 


(81) 


3z, 


’<♦9  '  Vi  "  Bi  W  * 
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(82) 


with  the  left  hand  sida  defined  as  sero  for  1  >  3 
and  for  energy  by 
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where  S1314  !•  th«  surface  cfitrfj  source. 

In  Eq.  79-63  the  symbol  I  represent#  e  general  dependent  variable, 
either  mass  or  energy  In  either  the  solid  or  gas  phase.  The  upper  and 
lower  signs  correspond  to  solid  and  gas  phaae*  respectively.  Volumetric 
and  msrface  sources  are  denoted  by  V  and  £,  respectively.  Other  details 
regarding  the  symbols  are  indicated  In  the  table  of  nomenclature. 


( 
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Appendix  A 

nrrERFACZ  cokdxttoks  for  equilibrium  ptroltsis 
or  on  or  mou  solid  phase  cokstitoemts 


Xhm  interface  species  uis  balances  for  pyrolyzlng  species  sre  provided 
by  Eq.  A 6-40  for  the  esse  trfiere  all  surface  processes  sre  controlled  by 
Arrhenius  type  rate  expressions.  At  the  other  estrone,  the  Interface  con¬ 
dition  nay  include  equilibrium  transfers  of  one  or  more  species  in  which  the 
rule  governs  the  relationship  between  concentration,  pressure,  end 
t**r*ratmre.  Consider,  for  example,  the  esse  of  equilibrium  transfer  of 
ingredisnt  10  (solid  fuel)  tdiile  Arrhenius  rate  expressions  govern  the  solid 
ox  id  liar  and  product  surface  conditions.  We  have  then  for  the  balance  of 
fuel: 


(A-I) 


which  may  be  rewritten  as 


(A-2) 


In  Eq.  A-2,  has  been  replaced  by  Im  ♦  mX09  thc  soemation  being  taken  over 
all  species  except  the  one  In  question  (fuel  in  this  case).  The  mass  frac¬ 
tion  Tj  (a  function  of  mole  fraction  and  molecular  weight  of  all  species) 
is  determined,  for  equilibrium  processes,  by  total  pressure  and  surface 
temperature  through  the  Clauiua-Clapeyron  relationship.  If  fuel  is  the 
only  component  undergoing  an  equilibria!  phase  change,  all  the  terms  in 
a  depend  on  rate  expressions  and  may  be  evaluated  froo  surface  tempera¬ 
ture  and  composition.  If  more  than  one  equilibrium  process  occurs,  a  re¬ 
lation  analogous  to  Eq.  A-2  applies  for  each  process,  leading  to  a  set  of 
simultaneous  boundary  conditions.  In  the  present  theoretical  analysis  and 
its  associated  computer  program,  equilibria:  processes  are  not  considered; 
the  condition  of  constant  surface  temperature  and  concentration  is  not 
obtainable  by  a  simplification  of  the  general  expressions  unless  convection 
and  surface  regression  are  neglected. 
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IMUim  AULTSIS  or  EFFEC1  Of  ABSOBFTIOfl  OF  EADIABT 
mXCT  BT  CAS  PHASE  SPECIES 


radiant  flux  la  that  tte  been  tuffara  no  ettcnnatlrn  In  passing  fro*  th« 
source  to  th«  surface.  Absorption  of  lnconlng  energy  by  gaseous  species 
has  tbs  dual  effect  of  adding  energy  to  the  gas  phase  while  partially 
blocking  it  free  the  propellant  surface.  Assuadng  that  only  the  evolved 
gases  contribute  to  energy  absorption,  that  energy  reflected  fron  the  sur¬ 
face  is  negligible,  and  that  equal  absorption  coefficients  apply,  the  gas 
phase  energy  equation  is  aodifiad  by  addition  of  the  ten 


In  Eq.  B-l,  L.  represents  the  unattenoated  source  flux  while  Y  is  the  sun 
of  all  s-eclss  produced  by  any  choalcal  or  physical  transfomation;  hence, 

r  -  2-T4  -  y9. 

He  nay  write 


(B-2) 


Also,  it  is  joted  that  all  new  species  nust  ultimately  originate  at  the 
surface,  so  that  at  tine  t 


(B-3) 


Carrying  out  the  differ  satiation  of  expression  (B-l)  leads  to  the  addition 
of  the  ten 


Prece&s  pge  fetal 
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to  Cq.  56.  Boro 


(B-5) 


Is  cho  flux  ot  the  surface  after  attenuation  by  the  entire  gas  coluon  and 
the  expression  (B4)  Is  the  energy  absorption  by  the  gas  at  position  x. 
Transformations  into  the  p  and  4  coordinate  systems  lead  to  the  additional 
terms 


(B-6) 


In  the  gas  phase  representation  of  Eq.  71  and 


in  Eq.  79. 

The  foregoing  Is  a  simplified  analysis  of  an  admittedly  complex  phenom¬ 
enon.  Expansion  of  the  ideas  to  include  other  effects  such  as  variable  6g 
and  gas  re-radlatlon  would  not  be  difficult  In  principle.  The  entire  matter 
la  usually  neglected  on  physical  grounds  and  in  addition  Is  omitted  in  the 
current  study  to  avoid  Introduction  of  more  complicated  numerical  analysis 
than  appears  warranted  by  probable  gains. 
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